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ABSTRACT 

Popular models for describing the luminosity-density profiles of dynam- 
ically hot stellar systems (e.g., Jaffe, Hernquist, Dehnen) were constructed 
with the desire to match the deprojected form of an i? 1 / 4 light-profile. Real 
galaxies, however, are now known to have a range of different light-profile 
shapes that scale with mass. Consequently, although highly useful, the above 
models have implicit limitations, and this is illustrated here through their ap- 
plication to a number of real galaxy density profiles. On the other hand, the 
analytical density profile given by Prugniel & Simien (1997) closely matches 
the deprojected form of Sersic R}/ n light-profiles — including deprojected 
exponential light profiles. It is thus applicable for describing bulges in spi- 
ral galaxies, dwarf elliptical galaxies, and both ordinary and giant elliptical 
galaxies. Here we provide simple equations, in terms of elementary and spe- 
cial functions, for the gravitational potential and force associated with this 
density profile. Furthermore, to match galaxies with partially depleted cores, 
and better explore the supermassive black hole / galaxy connection, we have 
added a power-law core to this density profile and derived similar expressions 
for the potential and force of this hybrid profile. Expressions for the mass and 
velocity dispersion, assuming isotropy, are also given. These spherical mod- 
els may also prove appropriate for describing the dark matter distribution in 
halos formed from ACDM cosmological simulations. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: kinematics and 
dynamics - galaxies: nuclei - galaxies: structure - stellar dynamics 

1 INTRODUCTION 

Both elliptical galaxies and the bulges of disk galaxies, hereafter collectively referred to 
as "bulges", possess a range of light-profile "shapes" that are well described by Sersic's 
(1963, 1968) R l ' n model (e.g., Caon, Capaccioli & D'Onofrio 1993; Young & Currie 1994; 
Graham et al. 1996; Graham 2001; Balcells et al. 2003). This model is a generalisation of de 
Vaucouleurs' (1948, 1959) i? 1//4 model which is known to be only appropriate for a subset of 
elliptical galaxies having Mb ~ —21 mag (e.g., Kormendy & Djorgovski 1989; Graham & 
Guzman 2003). The _R 1//4 model's limitation lies in the fact that it has only two parameters: 
a radial scale and a surface brightness scale. The actual curvature, or "shape" , of every R 1 ^ 
model is the same. 
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This restriction has carried over into computer simulations of "bulges". This is because 
the popular models that are used for describing the luminosity-density profiles of bulges, 
such as those from Jaffe (1983), Hernquist (1990), and Dehnen (1993, see also Tremaine 
et al. 1994), were created in order to reproduce an i? 1 / 4 light-profile when projected. All 
three of these density models have exactly the same outer profile slope, declining with radius 
as r~ 4 . Consequently, these models are limited in their ability to a) simulate the range of 
observed galaxy structures and b) quantify the evolution of these structures. 

Over the past two decades it has also become apparent that the most luminous (Mb < 
—20.5 mag) elliptical galaxies have partially-depleted stellar-cores (Kormendy 1985; Lauer 
1985; Ferrarese et al. 1994; Lauer et al. 1995; Gebhardt et al. 1996). On the other hand, 
the less luminous elliptical galaxies, sometimes referred to as "power-law" galaxies, have 
continuously curving Sersic profiles that continue all the way in to the resolution limit of 
HST images (Trujillo et al. 2004). A promising explanation, albeit not the only one proffered, 
for the depleted cores is that the giant galaxies formed from the dissipationless merger of 
two or more fainter elliptical galaxies. The subsequent gravitational ejection of stars by 
the inwardly spiralling supermassive black holes (SMBHs) — from the progenitor galaxies 
- scours out the core of the new galaxy (e.g., Begelman, Blandford & Rees 1980; Makino 
& Ebisuzaki 1996; Faber et al. 1997; Milosavljevic & Merritt 2001, 2005; Graham 2004). 
Furthermore, fundamental connections have been found between the SMBH mass and a 
bulge's: (i) magnitude (Magorrian et al. 1998; McLure & Dunlop 2002; Erwin, Graham 
& Caon 2002), (ii) light-profile shape (Graham et al. 2001, 2003a) and (in) kinematics 
(Ferrarese & Merritt 2000, Gebhardt et al. 2000). 

Given the above connections, it is obviously important to have models that are able to 
unite properly the domain of the black hole with the rest of the galaxy. In this paper we 
present a modification of Prugniel & Simien's (1997) density model that already matches the 
observed range of 'outer' profile shapes; our modification allows one to additionally model 
partially depleted cores. Moreover, and importantly, we also derive exact expressions for the 
potential and force in terms of (fast-to-compute) elementary functions, making it possible 
for simulations to explore the influence of cores and differing profile shapes. By setting the 
size of the partially depleted core to zero, the equations for the potential and force are 
applicable to the original Prugniel-Simien density profile, for which no previous expressions 
existed. 

In subsection 2.1 we introduce the density model, complete with power-law core, while 
subsection 2.2 provides the equations for the potential and force. Expressions for the pro- 
jection of the density model can be found in subsection 2.3. In section 3 we fit a number 
of popular density models to the luminosity-density profiles of real galaxies and compare 
the results with the fit from the Prugniel-Simien density model and our modification of this 
model. A discussion of relevant issues is given in section 4, and section 5 provides a summary 
of the main results. Appendix A provides the derivation for various expressions pertaining 
to the velocity structure of the density profiles, including the circular velocity, spatial and 
line-of-sight velocity dispersion. Appendix B provides the derivation of the equations given 
in section 2. 

2 THE MODEL 

In this section we introduce a family of expressions associated with the spatial density 
profiles of spherical stellar systems having Sersic-like profiles with optional power-law cores. 
We derive exact expressions for the potential and force and outline how the corresponding 
surface brightness can be computed numerically (see FigJTJ). With computer processing speed 
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Figure 1. Flow chart for the density model: the projection of the analytically prescribed density is computed numerically, 
while the potential and force are obtained exactly. BT 2-22 refers to Binney & Tremaine (1987, their equation 2-22). 



in mind, in addition to the pressing need for more flexible galaxy models, the equations 
consist only of analytical terms and simple elementary functions. 

For the sake of parameter identification, we first introduce Sersic's (1963, 1968) model 
for describing the projected, radial intensity profiles of galaxies, such that 

I(R) = I e- b W R ^ 1/n , (1) 

where Iq is the central intensity and R c is the (projected) effective half-light radius. The 
parameter n describes the curvature of the profile (see, e.g., Ciotti 1991). The term b is not a 
parameter, but instead a function of n and chosen to ensure R e contains half the (projected) 
galaxy light. It is obtained by solving the equation r(2n) = 2 7(272, b), where 



T(a) = l°° e'H^dt and 7(0, x) = / VV^Mt, a>0 (2) 
io Jo 

are the complete and incomplete gamma functions, respectively. Although we have chosen to 
numerically compute the exact value for b because fast codes are available for the computa- 
tion of gamma functions, a good approximation of b for 0.5 <n< 10 is 2n — 1/3 + 0.009876/ra 
(Prugniel & Simien 1997, see MacArthur, Courteau, & Holtzman 2003 for smaller values of 
n). Using the substitution x = b(R/ ' R c ) l l n , the total luminosity from equation (JTJ) is given 
as 

POO 

Ltot = / I(R')2irR'dR' = 27iI R 2 c nb- 2n T(2n). (3) 
Jo 

For a sample of 250 dwarf elliptical and ordinary elliptical galaxies spanning —13 > 
Mb > —23 mag, Graham & Guzman (2003) have shown how profile shape and surface 
brightness vary with galaxy magnitude. They found Mb = — 9.41og(ri) — 14.3, and Mb = 
(2/3)/io — 29.5 mag — until the presence of cores in galaxies brighter than ~-20.5 I?-mag. 
The latter phenomenon makes it more appropriate to use /i e , the surface brightness at R e , 
rather than the central surface brightness fiQ = — 2.51og(/o)- Doing so gives the expression 
M B = (2/3)(/i c — 1.0866) — 29.5 mag. Typical values are given Tabled] 

A compendium of expressions related to Sersic's R l ^ n model can be found in Graham & 
Driver (2005). 
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Table 1. Typical galaxy masses and central surface brightnesses /jo = — 2.5 log(7b) (using the inward extrapolation of the outer 
Sersic profile in the case of "core" galaxies) associated with a range of Sersic indices n. The mass estimates have come from the 
absolute B-band magnitudes in Graham & Guzman (2003) using M b q=5A7 (Cox 2000) and M/L s =5.31 (Worthey 1994, 
assuming a 12 Gyr population with Fc/H=0). The numbers are only indicative, with intrinsic variance playing a role amongst 
real galaxies. 
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2.1 Density 

Generalising an expression from Mellier & Mathez (1987) that approximated the spatial, 
i.e. not projected, density profile of the _R 1//4 model, Prugniel & Simien (1997) provide 
an analytical approximation to the density profile of the R}l n model. Lima Neto, Gerbal 
& Marquez (1999) showed that this spherical model is accurate to better than 5% over 
the radial range 10 _2 -10 3 i? c . An even more accurate approximation for both spherical and 
triaxial stellar systems with R 1 ^ light-profiles was developed by Trujillo et al. (2002), with 
an accuracy better than 0.1%. However, although analytical, the latter expression is not 
particularly simple. Impressively, an exact solution to the deprojection of the R}l n model, 
which obviously includes the R 1 ^ case, was given in terms of Meijer G functions by Mazure & 
Capelato (2002); but again these equations are somewhat complicated 1 . Moreover, although 
these models can properly treat the range of outer profile shapes observed in real galaxies, 
they can not additionally allow for the presence of partially depleted cores. In fact, we are 
unaware of any density model capable of simultaneously describing the range of structure 
observed in both the inner and outer regime of galaxies with "cores" . 

The 3-parameter (p , -Re, n) density profile of Prugniel & Simien (1997; their equation B6) 
is relatively simple and can be written as 



(4) 



P0 L 2R c T(n{3-p)) ' 

where r is the spatial radius and po is the normalisation such that the total mass from 
equation (J3J see Appendix A) equals that from equation (JTJ) . We adopt Lima Neto et al.'s 
(1999) estimate, or rather the updated value given in Marquez et al. (2000), for the term 
p, for which a high-quality match between the exact, deprojected Sersic profiles (solved 
numerically) and the above expression is obtained when p — 1.0 — 0.6097/n + 0.05563/n 2 , 
for 0.6 < n < 10 and 10~ 2 < R/R e < 10 3 . The quantity M/L = T is the mass-to-luminosity 
ratio, which is typically taken not to depend on galaxy radius but may of course do so. The 
density at r = R e is simply poe~ b . Figure (0) shows the behaviour of this density profile for 
different values of profile shape n. Expressions and figures for the (enclosed) mass profile, the 
circular velocity, and the spatial and line-of-sight velocity dispersion are given in Appendix 
A for the spherical case. 

Adding an inner power-law with slope 7, which is not to be confused with the incomplete 
gamma function 7(0, x), we obtain the new model 



Exact analytical expressions for the mass, gravitational potential, total energy and the central velocity dispersion are also 
presented in Mazure &; Capelato (2002). Exact numerical expressions arc given in Ciotti (1991). 




Figure 2. Core-less galaxies. Panel a) projected intensity / (exact numerical solution); panel c) potential (equation IS1 divided 
by <E>(0.01-fic)); panel d) force (equation 1141 : and panel e) pseudo potential (equation 1151 divided by |<J>(0.01ilc)| with h 2 = 
10 — 6 |4>(0.01-f? c )|) associated with the Prugniel-Simien density profiles (equation |1J shown in panel b) for varying values of the 
profile shape n: n = 0.5 (solid lines), n = 1 (dotted), n = 2 (dashed), n = 4 (dash-dot), n = 10 (double— dash). The radius and 
density are normalised such that R c = 1 and p(0.01i?, c ) = 1. 



p{r) =f/[l+ (^) a y /a {[(r a + r b a )/R c a r P/a e- b i {ra+r Z )/R °y /na 
p> = Pb 2&-r>/« (^) P e b ( 2l/Q ^) 1/n . 



(5) 



The break radius, r b , denotes the transition where the profile changes from one regime to 
the other, with p b the density at this radius. The parameter a controls the sharpness of 
the transition. For r << r& equation (jSJ) tends to a power-law with slope 7. For r >> r& 
equation © reduces to equation (jlj). 

Modelling the light-profiles of luminous elliptical galaxies, Trujillo et al. (2004) have 
shown that the transition from the inner "core" to the outer Sersic profile is probably rather 
sharp. Motivated by this result, we have chosen to make the transition between the inner 
power-law and the outer profile sharp by considering the a — > 00 limit. This reduces the 
above expression to a 5-parameter model capable of describing the entire radial extent of 
spherical stellar systems with power-law "cores" , and can be written as 



P(r) = Pb 



(a) 7 h (r h -r) + p (£)~ P e~ b ^ R ^ 1/n h(r - r h ) 



(6) 



where h(x) is the Heaviside step function such that h—1 if x > and h=0 if x < 0. 

These spherical models are illustrated in Fig|2]e-h). In each panel a range of values of 
n has been used, while the value of 7 increases sequentially from to 1, 1.5 and finally 2. 
Depending on the parameter combination, one may have either a partially depleted core 
or a central excess — possibly representative of a nuclear star cluster or adiabatic growth 
around a supermassive black hole (e.g., Merritt 2004, and references therein). The break 
radius has been set to 0.01i? e , which a) matches the observed core-radii values of 0.01 to 
0.02i? c found by Trujillo et al. (2004, their table 2), and b) only uses the density profile 
from Prugniel & Simien (1997) over the radial range where it provides a high-quality match 
to a deprojected Sersic profile. For reference, the nuclear star cluster "half-width at half- 
maximum" (HW H M) values divided by their host bulge R e values in nucleated dwarf 
elliptical galaxies is 0.02 to 0.04 (Graham & Guzman 2003, their table 2). 



2.2 Gravitational Potential and Force 

We derive an expression for the potential (assuming spherical symmetry) using the inverted 
form of the Poisson equation (e.g., Binney & Tremaine 1987, their equation 2-22): 
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Figure 3. Core galaxies. Panels a-d) show the projected intensity / (equation 1181 exact numerical solution); panels i— 1) show the 
potential (equation [S]dividcd by <I>(r = r^)); panels m-1) show the force (eouation ll4l ; and panels q— t) show the pseudo potential 
(equation I15l dividcd by |<J>(r = rjj)| with h' 2 = 10~ 6 |<E>(r = r^)\) associated with the density profiles (equation 151 in panels e— h), 
for varying values of the outer profile shape n: n = 0.5 (solid lines), n = 1 (dotted), n = 2 (dashed), n = 4 (dash-dot), n = 10 
(double-dash); and varying central cusp slope 7. In these figures the scale radius R e has a value of 1, as does the scale density p;, 
which occurs at a radius = 0.01i? c . 



$(r) = -47rG 



- J p(r)r 2 dr + j 
r 



After some algebraic manipulation (see Appendix B), one obtains 

- -4vrG { \ Jl(r) + ¥ r) + Ll{rb) if r " r? " 
where 

Ji(r) = p 6 r 6 V-7(3 - 7) if 7 <3 

jL. ( rfo 2-7 _ r 2-7) if 7 ^ 2 

if 7 = 2 ' 



L x (r) 
and 



p b pR e 2 nb n(p - 2) T (n(2 -p),b (jffj 
Pb pR e 3 nb n{p - 3) 



POO 

r(a,x) = / e-H a - l dt, a > 



(7) 

(8) 

(9) 
(10) 

(11) 
(12) 

(13) 
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is the complement to the incomplete gamma function shown in equation (J2J). Although J x 
diverges when 7 = 3, this is unlikely to be a problem because such steep inner profile 
slopes are not observed in real galaxies (e.g. Gebhardt et al. 1996; Milosavljevic et al. 2002; 
Ravindranath, Ho, & Filippenko 2002). FigOHHl show the gravitational potential of this 
model as a function of spatial radius r. The case when the break radius equals zero is shown 
in Fig. (j2t). 

The central potential $(0) = -47rG[J 2 (0) + £i(r 6 )], for 7^2. When r b = 0, $(0) = 
— 4-7rGLi(0), with pf,p = p - 

The corresponding radial force is computed by differentiation with respect to the spatial 
radial coordinate r. After some algebraic cancellations (see Appendix B), the force is found 
to be 



We remind readers of the logarithmic radial scale used in Fig. |3J and caution that the 
gradient to the curves shown in Fig|3HIIl should be interpreted with care when comparing 
them with the force shown in Fig|Hfri-|np- The case when the break radius equals zero is 
shown in Fig. (j2Ji). 

In order to carry out orbital integration, one needs to compute the force along an orbit. 
For a fixed break radius r^, the term in the first half of the expression L 2 (r) is constant and 
therefore only needs to be computed once at the beginning of a simulation. Therefore, the 
incomplete gamma function T(a, x) need only be called once per evaluation of the force (the 
term in the second half of L2 (?"))> anc ^ no ^ a ^ a ^ ^ the star is inside the break radius, i.e. 
r < ri,. 

The abrupt change in the force — whose nature is not only dependent on the value of the 
inner profile slope 7 and the Sersic index n but also on the ratio r&/ R c — is not a problem 
for constructing models because the force remains both negative and continuous at r&. Orbit 
stability can be checked using the "pseudo" potential ^(r), sometimes also referred to as 
the "effective" potential (Landau & Lifshitz 1976; Goldstein 1980), and given by 



where h is the angular momentum per unit mass. Uisng h 2 = 10- 6 |$(0.01i? e )|, and normal- 
ising \l/(r) by dividing by |$(0.01i? e ) |, in order to highlight orbits near the radius 0.01i? c , 
FigEh-S reveals that the abrupt change in density and force at r& does not result in a 
double minimum for the pseudo potential — which would be indicative of an unstable con- 
figuration. Instead, for a given orbital energy less than zero, particles remain bound within 
the curve defined by the pseudo potential. 

In summary, we have a set of equations consisting of analytical expressions, and one 
special function for which fast computer codes exist, that can a) simulate the range of 
observed light-profile shapes and b) model partially depleted cores and/or certain additional 
nuclear components. 

2.3 Projected Surface Brightness Profile 

The projection of the above density profile, i.e. the intensity-profile or "light-profile", is 
computed by solving the Abel integral 




(14) 



M/( r ) = $( r ) + 



0- 

2^2' 



(15) 




(16) 
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where R is the projected radius and T is the mass-to-luminosity ratio. If equation (J^l is 
rewritten as p(r) = pi(r) + p 2 (r), where 

Pi(r) = Pb(-) h(r b -r), and 



r 

p 2 (r) = p b pi^y\-Kr/^ n h{ r_ r ^ (17) 
then the corresponding intensity can be expressed as 

\ I ^Pwdr + 4 / ^p^dr = h(R) + I 2 (R- r b ) if R < r b , 
I(R) = { (18) 
if R j&dr = I 2 (R) if R>r b . 

As shown in Appendix B, the first term Ii(R) is such that 



r\ -R 2 if 7 = 0, 

(19) 



In rb+ Vf^ if7 = l 



h(R) = 

T R~ x sin" 1 - # 2 /r 2 if 7 = 2, 

^R^Bi^tf/r^ Q, otherwise, 

with i? the incomplete Beta function defined as 

X 

B x (y,z) = [ u^\l-uf~ l du. (20) 



The only expressions we do not provide an exact equation for are hiR) and I 2 (R;r b ). 
However, because Prugniel & Simien (1997) devised equation (jU) to match the deprojected 
form of the Sersic light-profile, it makes sense to use the R l ' n model as a suitable approxi- 
mation. One therefore has that 

I 2 {R) « I e e b e- b W R ^ /n (21) 

with I c the (projected) intensity at the (projected) radius R e . The value of I c can be set in 
terms of the density model parameters if one applies the condition that the total mass given 
by equation (|2*T1 see equation EJ) equals the total mass from equation (jU see Appendix A), 
yielding 

2e- b R ePb pT(n(3-p)) 

Tr(2n)H 1 -f) ' { } 

with p given in equation (JHJ). Over the radial interval 10~ 2 < R/R c < 10 2 , the maximum 
difference in surface brightness between this approximation for I 2 (R) and the exact value 
is about 0.1 mag arcsec -2 if n > 2, and only 0.04 mag arcsec -2 when n = 4 (Fig. HJ). For 
a Gaussian profile (n = 0.5) the agreement is very good, while for an exponential profile 
(n = 1) the match at large radii is rather poor — increasing to ~0.2 mag arcsec -2 at 
R/R c = 100. 

The only term requiring numerical evaluation is I 2 (R; r b ), although, if desired, one may 
choose to additionally evaluate I 2 {R) numerically rather than using the analytical approx- 
imation given in equation (J21j) . Given that the intensity profile is likely only computed at 
the end, or at a few intermediate stages, of a simulation, the computation time involved in 
deriving the light-profile is not a concern. 

The exact (numerically computed) projected light-profiles associated with the density 
models in Fig. (J2b) are shown in Fig. (j2t), and those in Fig|3kH2h are shown in FigO^HSti- 
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Figure 4. The relative error between the term l2(R) given in equation 1181 and the approximation given in equation 1211 is 
shown as a function of radius, normalised by the effective radius R e , for varying light-profile shapes n. 

In passing we note that when 7 = and R < r&, the slope of the surface brightness profile 
is close to but not exactly zero (FigJHk,). We also note that the value of I e and R e are stable 
against increases in and changes in 7 because they pertain to the outer, undisturbed profile 
when there is a partially depleted core, and to the underlying host galaxy profile when a 
central flux excess is present. As such, when r& > 0, I e and R c are not the total galaxy half- 
light values. In practice, the stellar flux deficit is only ~0.1% for "core" galaxies (Graham 
2004) and the stellar flux excess ~1% in nucleated dwarf elliptical galaxies (Graham & 
Guzman 2003), and so the actual discrepancy is not greater than 

3 APPLICATION TO REAL GALAXY PROFILES 

Over a quarter of a century ago, de Vaucouleurs & Capaccioli (1979) demonstrated that, 
when the inner 10 arcseconds were excluded, the R 1 ^ model provided a very good fit to the 
light-profile of NGC 3379 (Mb ~ —20 mag). Since then, Capaccioli and his collaborators 
have shown that different galaxies can be equally well fit, but only by using models with a 
very different curvature (e.g., Caon et al. 1993). Bertin, Ciotti & Del Principe (2002) have 
shown R l l n fits to the circularised light-profiles of four (illustrative) galaxies from Caon 
et al. (1993). Here we deproject these galaxies' major-axis light-profiles and fit the Jaffe, 
Hernquist and Dehnen models, and explore how well they describe the luminosity-density 
profile in comparison with the model of Prugniel & Simien (1997) given in equation 
In addition, we include the i/5T-resolved dwarf elliptical galaxy 2 LGC 47 (Stiavelli et al. 
2001), and the "core" galaxy NGC 3348 (Rest et al. 2001; Trujillo et al. 2004) to which we 
apply our new core-density model (equation 

3.1 Sersic and core— Sersic model 

FigEl shows the Sersic (1963, 1968) model applied to the ground-based, major-axis, B- 
band light-profiles of NGC 1379, 4458, 4374, and 4552, and applied to the HST-b&sed, 
major-axis, J-band light-profile of LGC 47. The data have come from Caon et al. (1993) 
and Stiavelli et al. (2001), respectively. Also shown is the core-Sersic model (Graham et al. 
2003b; Trujillo et al. 2004) applied to the i/ST-based, major-axis, i?-band light-profile of 



2 LGC: Leo Group Catalog. 
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Figure 5. Best-fitting Sersic model to the major-axis, S-band light-profiles of NGC 1379, 4458, 4374 and 4552 (data 
from Caon et al. 1990; Caon et al. 1993), and the major-axis, /-band light-profile of LGC 47 (data from Stiavelli et al. 
2001). The best-fitting core— Sersic model to the major-axis, _R-band light-profile of the "core" galaxy NGC 3348 is also 
shown (data from Trujillo et al. 2004); the inward extrapolation of the outer Sersic profile is shown by the dashed curve. 
The model parameters are given in Table The effective half-light radii are marked with an arrow. 
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NGC 3348. Table El shows the best-fitting parameters. For the first four galaxies, the values 
agree with those reported in Caon et al. (1993), with the exception that we find a slightly 
smaller value of n for NGC 4552. However, the difference between an n = 12 and n = 14 
profile is minimal. For LGC 47, Stiavelli et al. (2001) presented fits to the geometric mean 
(r = \^ab), V-band light-profile, finding n ~ 1.5. Due to possible colour gradients, and 
ellipticity gradients, it is expected that our fit to the major-axis, /-band light-profile may 
be slightly different: nonetheless, we derived a similarly small value of n = 1.1. NGC 1379 
and NGC 4458 have values of n = 2.0 and 2.6 respectively. NGC 4374 and NGC 4552 
have values of ~ 8 and ~ 12, considerably greater than 4. On the other hand, NGC 3348 
has a value of n close to 4, but possesses a distinct core. Although it is likely the large 
galaxies NGC 4374 and NGC 4552 may also possess a partially depleted core, the inner few 
arcseconds have been excluded from their profile due the effects of seeing. 



3.2 Density profiles 

In order to obtain each galaxy's deprojected light-profile, a non-parametric deprojection 
was used when solving the appropriate Abel integral. The best-fitting light-profile models 
were however used to extrapolate the observed light-profiles to infinity. This simply meant 
using the best-fitting Sersic model to extend the observed data to larger radii in order for 
one to compute the deprojected light-profile, i.e. the density profile. The actual choice of 
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extrapolation only affects the outer luminosity-density profile slightly. Due to the coarser 
radial sampling of the ground-based profiles than the HST profiles, their deprojected profiles 
are somewhat noisier. We applied box-cart smoothing with four different box sizes and 
refitted the various density models; in every instance the best-fitting model parameters 
where practically identical, suggesting such noise is not an issue. 

In the case of the dwarf elliptical galaxy LGC 47, its observed (exponential) light-profile 
only went out to 1 R e , not enough to properly show how such stellar distributions decline at 
large radii. For this galaxy only, we model the density profile to radii larger than observed. 
Given that galaxies, and bulges, with exponential profiles do exist, we felt that the analysis 
of such a profile would be of value. 

The luminosity-density profiles (FigEJ^J) have been calibrated in units of Lq pc -3 
using the distance moduli provided in Tonry et al. (2001) for the Caon et al. sample, and 
a distance of 41 Mpc for NGC 3348 and 11.5 Mpc for LGC 47, taken from their respective 
papers. We used solar absolute magnitudes of M B =5.47, M R =4.28 and M/=3.94 (Cox 2000). 
A Hubble constant H = 75 km s -1 Mpc^ 1 was used. 

The value of A associated with the residual profiles shown in FigloTfTUl is not quite the 
rms error, but is computed using the expression 



where m is the number of data points, 5i is the i th residual and k is the number of parameters 
in the fitted model. Because the different models have different numbers of free parameters, 
this provides a more appropriate measure for model comparison. 

3.3 Jaffe model 

Fig presents the best-fitting (2-parameter) Jaffe (1980) model, p(r) = [4p(a)](r/a) _2 (l + 
r/a)~ 2 , to each galaxy's luminosity-density profile. The scale-length is denoted by a, and 
p(a) is simply the density at r = a. The position of the transition radius, i.e. the scale-length, 
where this model changes from an inner power-law slope of —2 to an outer power-law slope 
of —4 is marked with an arrow in FigEl in the case of LGC 47 and NGC 3348, the Jaffe 
model is clearly inadequate to describe the stellar profile. We also note that the best- 
fitting models have transition radii well beyond the observed radial extent of the galaxies. 
The hump-shaped residual profiles (upper panels) for the low n galaxies NGC 1379 and 
NGC 4458, and the bowl-shaped residual profiles for the high n galaxies NGC 4374 and 
NGC 4552 are clear indications that the Jaffe model is unable to match the global curvature 
in these galaxies' stellar distributions. 

3.4 Hernquist model 

In FigJ7|one can see the best-fitting (2-parameter) Hernquist (1990) models, p{r) = [Sp(a)](r / a) 
r/a)~ 3 . This model is also incapable of describing the core galaxy NGC 3348, and the 
hump/bowl shaped residual profiles of the other galaxies largely reflects the failure seen 
with the Jaffe model. We do however note that in the case of NGC 4458, the (non-core) 
galaxy with the closest Sersic index to a value of 4, the fit is quite acceptable. 

3.5 Dehnen model 



The (3-parameter) Dehnen (1993) model, p(r) = [2 4 7 p(a)](r / a) 7 (1 +r/a) 7 4 , which sub- 
sumes the first two models as special cases, does noticeably better. This is due to the fact 
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Figure 6. Best— fitting Jaffo model to the deprojected light-profiles of the galaxies in Fig|F] The luminosity-density profile 
u(r) = p(r)/(M/L) is in units of Lq b pc" 3 for NGC 1379, 4458, 4374 and 4552, L & Ji pc~ 3 for NGC 3348, and L Q l 
pc — 3 for LGC 47. The model parameters are given in Tablc|2] The arrow denotes the size of the scale radius a. 
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Figure 7. Same as Fig|H] except that the best-fitting Hernquist models are shown. 



that it has an additional (third) parameter over the previous two models that actually 
measures the inner power-law slope, rather than assuming some fixed value. Obviously the 
Dehnen model fails to fit the detailed profile of NGC 3348. In a sense, it also fails to match 
NGC 4552, setting the transition radius to 500" — an arbitrary large upper limit in our 
code — well outside of the observed radial range. The reason for this large transition radius 
is however understood. Graham & Driver (2005) show that for large values of n, the Sersic 
model tends to a power-law such that I(R) ~ R~ 2 ] which results in p(r) ~ r~ 3 . NGC 4552 
has the largest value of n in our sample, and the optimal Dehnen model is simply a (single) 
power-law with slope 7. 

In the case of the lower n galaxies (LGC 47 and NGC 1379), one can see that the 
Dehnen model — constructed to match i? 1 / 4 light-profiles — does not decline as quickly as 
the observed structures. Due to the fact that we used an exponential-like model (n—1.1) 
to extrapolate the light-profile of LGC 47, this eliminates the possibility that the observed 
mismatch at large radii may be unique to this galaxy. Rather, the mismatch reflects the 
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Table 3. Best-fitting parameters from a range of density models. The units are Lq pc 3 for p and arcseconds for a. 
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well known fact that exponential light-profiles have outer density profiles that decline more 
quickly with radius than r~ 4 . Not surprisingly, an almost identical result to that seen in 
Fig- El for LGC 47 is obtained when we deproject a pure n — 1 exponential light-profile. As 
a rule, the lower the value of the Sersic index n, the steeper the outer decline in density and 
the greater the mismatch with the Dehnen model. 

3.6 Prugniel— Simien density profile, and our adaptation for galaxies with 
power— law cores 

FigElpresents the best-fitting (3-parameter) model of Prugniel & Simien (1997). Aside from 
NGC 3348, the fits are obviously rather good, with only NGC 4552 suggesting the presence 
of additional fine structure that has been missed. In fact, Caon et al. (1993) classified this 
as a lenticular galaxy; the other galaxies are elliptical, i.e. they do not have an embedded 
large-scale disk. With the exception of the core galaxy NGC 3348, every fit is equal to or 
better than those obtained with the (3-parameter) Dehnen model. This is especially the 
case for the low n galaxies LGC 47 and NGC 1379. 

FigEl presents the new (5-parameter) model (equation EJ) applied to the core galaxy 
NGC 3348. The model performs well and is the only density model we are aware of that 
can fit galaxies possessing partially depleted cores. The best-fitting parameters are given in 
Table El Although the luminosity-density at a spatial radius r equal to R c is not a formal 
parameter of the new model, it is equal to 10~ L26 Lq r pc -3 for NGC 3348. 

We have also included two additional "core" galaxies: NGC 2986 and NGC 4291. The 
best-fitting core-Sersic parameters derived here for NGC 4291 are the same as those given in 
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Figure 9. Same as Fig|H] except that the best-fitting model of Prugniel & Simien (1997), as given by equation El . is 
shown here. The arrow marks the radius R c . 
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Figure 10. New density model (equation|^J applied to the core galaxies NGC 3348, 2986 and 4291. The arrow marks the 
model radius R e . 



Trujillo et al. (2004), as was the case with NGC 3348. However, our core-Sersic parameters 
for NGC 2986 are somewhat different; the reason is due to the fact that we exclude the 
(highly deviant) outermost two data points from the light-profile given in Trujillo et al., 
which may be due to poor sky-subtraction. Although the core-parameters did not change 
much, the value of n increased by 25% and R e increased from 44 to 76 arcseconds. 

NGC 2986 and NGC 4291 were not shown in the previous figures because they had 
residual profiles very similar to that seen for NGC 3348. They are included here so one can 
see the new model's applicability to other core galaxies. Interestingly, the residual profiles 
seem to suggest that the core density profiles have only approximately a power-law structure. 



4 DISCUSSION 

Just like the R 1 ^ 4 model, the R x l n model was offered as an empirical fitting function. A 
physical basis for its origin is not widely recognised. Marquez et al. (2000, 2001, and ref- 
erences therein) have argued it is a natural result from the (near) conservation of specific 
entropy. Other theoretically motivated models, based on the statistical mechanics of par- 
tially complete violent relaxation, also lead to density profiles that match deprojected Sersic 
R x l n profiles, at least for 2.5 < n < 8.5 (Trenti & Bertin 2005, see also Hjorth & Mad- 
sen 1995). We do not, however, attempt to solve this question here, but do note that the 
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Jaffe, Hernquist and Dehnen models, like the new model presented here, are simply useful 
empirical fitting functions. 

FiglHl and [7| clearly reveal that the Jaffe and Hernquist models fail to describe the stellar 
distribution of low- and high-n galaxies, relative to n — 4. This is not particularly surprising 
as they were only designed to match galaxies with an n = 4 light-profile. The Dehnen model 
does better, but fails to model the lower n galaxies in our sample. This can be seen in the 
outer part of the profile for NGC 1379 and LGC 47 in Fig- El where the density declines more 
rapidly with radius than r -4 . At the high-n end, the (projected) intensity profile I(r) — > r~ 2 
as the Sersic index n — > oo (Graham & Driver 2005), and the density profile tends towards 
a power-law with p(r) — > r -3 . For reference, the isothermal model has p(r) oc r -2 . Ignoring 
the presence of partially depleted cores, for galaxies with values of n greater than ~8, a 
single power-law can approximate their density profiles and the Dehnen model appears to 
work well, albeit by setting the scale radius, a, to very large radii. 

Zhao (1996) developed a 5-parameter generalisation of the Dehnen model such that 
both the inner and outer power-law slope can be adjusted, along with the radius, density 
and sharpness of the transition region. The model was first introduced by Hernquist (1990, 
his equation 43) and has the same structural form as the "Nuker" model (e.g., Lauer et 
al. 1995). However, as explained in Graham et al. (2002, 2003b), such double power-law 
models are not appropriate for describing profiles with obvious logarithmic curvature, that 
is, profiles without inner and outer power-laws but whose slopes continuously vary as a 
function of radius — which is the case for the luminosity-density profiles of most elliptical 
galaxies and bulges. Although in the absence of a partially depleted stellar core this model 
should be able to approximate real luminosity-density profiles, the parameters themselves 
are highly sensitive to the fitted radial range. This is because the parameters simply adjust 
themselves in order to match the curvature in whatever part of the profile one includes 
in the fit, providing inner and outer power-law slopes that are a product of the radial 
range sampled. Consequently, this model is not explored here. The 3-parameter model of 
Prugniel & Simien (1997) was, however, designed to describe profiles with curvature and 
does therefore not suffer from such a problem. 

The density profile of Prugniel & Simien (equation |4} may prove helpful for describing 
gravitational lenses (e.g., Cardone 2004; Kawano et al. 2004) and the dark matter distribu- 
tions of halos built in ACDM cosmological simulations. Demarco et al. (2003) have already 
observed the hot gas in galaxy clusters to have a (projected) Sersic distribution. Intriguingly, 
the curved nature of simulated dark matter profiles has resulted in the recent discovery that 
Sersic's model describes them better than the generalised NFW 3 double power-law model 
with inner and outer slopes of —1 and —3, respectively (Merritt et al. 2005). Whether or 
not they were aware of the fact, Navarro et al. (2004; their equation 5) and Cardone, Piedi- 
palumbo & Tortora (2005; their equation 10) presented Sersic's model as a leading candidate 
to describe the density profiles of dark matter halos. It would be of interest to investigate 
whether the generalisation of the Mellier-Mathez (1987) model presented by Prugniel & 
Simien (1997) (equation 0] in this paper) may be useful, if not even more appropriate, for 
such studies. 

The new density-potential pair should also prove fruitful for simulations of bulges in 
disk galaxies. To date, such studies have usually been performed using models which are 
inappropriate for describing the majority of bulges. This situation has arisen from improved 



3 If cosmological simulations of dark matter halos indeed have density profiles which are logarithmically curved, i.e. if they 
have slopes which vary continuously with radius (e.g. Navarro et al. 2004; Hansen & Moore 2005), then the issue of whether 
the central cusp slope is -1 or -1.5 or smaller may need to be reconsidered. 
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observations which have shown bulges typically have Sersic indices n ranging from 0.5 to 3, 
rather than 4 (e.g., Andredakis & Sanders 1994; de Jong 1996; Graham 2001; MacArthur et 
al. 2003). 

It is hoped that the models presented here will additionally enable one to perform a num- 
ber of simulations testing various observational results and theoretical problems associated 
with the influence of SMBHs on central density cusps and the connection with the global 
stellar structure. The influence of SMBHs on the stellar dynamics (e.g., Baes, Dejonghe & 
Buyle 2005) can also be explored as a function of various physical parameters. One way this 
can be achieved is by investigating the restrictions that self-consistency requirements pose 
on galaxy density profiles (e.g., Terzic 2003). Other numerical simulations may investigate 
the effects of binary SMBHs at the centers of galaxies (e.g., Graham 2004). Gravitational 
scattering (e.g., Makino & Funato 2004), as well as resonant chaotic phase mixing between 
the frequencies of the binary orbit and the natural frequencies of stellar orbits, can cause a 
significant redistribution of mass extending well beyond the radius of the binary (Kandrup 
et al. 2003; see also Merritt & Poon 2004). One could also address the similarities and dif- 
ferences between the nature and efficiency of chaotic phase mixing in the new model (both 
time-independent and subjected to time-periodic pulsations) and previous studies that have 
used somewhat more limited potentials (e.g., Kandrup & Mahon 1994, Merritt & Valluri 
1996; Siopis & Kandrup 2000, Terzic & Kandrup 2004). In addition, one should be able 
to extend past work which has explored the evolution of galaxies using models that had a 
restricted range of outer galaxy structure (e.g., Merritt & Fridman 1996; Poon & Merritt 
2002, 2004). Furthermore, it is hoped one will explore new connections between the outer 
galaxy structure, i.e. profile shape, and the properties of the core, such as the inner cusp 
slope, core size, and black hole mass. 

5 SUMMARY 

Our application of the (2-parameter) Jaffe and Hernquist models to a sample of luminosity- 
density profiles taken from real galaxies reveals that these models are inadequate for de- 
scribing the observed range of stellar distributions. The (3-parameter) Dehnen model does 
better, providing a good match to the density profiles of large early-type galaxies with Sersic 
indices around 4 and greater, but it fails to adequately match the deprojected light-profiles 
of galaxies with low Sersic indices, which would encompass dwarf elliptical galaxies and most 
bulges in disk galaxies (e.g., Balcells et al. 2003, and references therein). If one wishes to 
explore the hierarchical merging theory, in particular the potential build up of large elliptical 
galaxies through the collision of lesser elliptical galaxies, then none of the above models are 
good starting points. 

These failures are not particularly surprising because all three of these models were 
constructed to have outer density profiles that fall off as r~ 4 — in order to match the decline 
in density of a deprojected i? 1 / 4 profile. As such, they are highly useful but nonetheless 
limited models. The (3-parameter) density profile of Prugniel & Simien (1997), however, 
does not have the above constraint, and consequently does much better, accurately matching 
the density profiles of luminous elliptical galaxies with n > 4 and faint ellipticals with n < 4, 
including exponential (n = 1) profiles. This model may even prove applicable to the density 
profiles of ACDM-generated dark matter halos. 

We have derived exact expressions for the gravitational potential and force associated 
with the latter density profile. The use of this family of expressions for simulations of el- 
liptical galaxies and bulges in general, and explorations of how their structure may evolve 
under various circumstances, should enable projects of a nature previously prohibited by 
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the former class of models. Moreover, the equations for the potential and force contain only 
one elementary function and are otherwise analytical; they are therefore fast to compute. 

We have also developed what we believe to be the only density profile capable of simulta- 
neously matching both the nuclear and global stellar distribution in galaxies having partially 
depleted cores. Application of this (5-parameter) profile to NGC 3348 reveals a good fit (rms 
~ 0.03 dex) over 3 orders of magnitude in radial range and 6 in luminosity density. Similar 
results are obtained from the core galaxies NGC 4291 and NGC 2986. Furthermore, the 
associated equations for the potential and force are also derived here using only the incom- 
plete gamma function and analytical terms. Expressions to derive the enclosed mass, the 
spatial and projected velocity dispersion, and the projected intensity profile are additionally 
provided. 

6 ACKNOWLEDGMENTS 

We are happy to thank Chris Hunter, Andi Burkert and Agris Kalnajs for their helpful 
comments and suggestions. We also wish to thank the referee, Eric Emsellem, for his com- 
ments on the text. We are grateful to Nicola Caon for providing us with the light-profiles 
for NGC 1379, 4374, 4458 & 4552, to Peter Erwin for kindly providing us with the light- 
profile for NGC 3348 and to Massimo Stiavelli for supplying us with the light-profile for 
LGC 47. This research was supported in part by NSF grant AST-0307351, Department of 
Energy grant G1A62056 and by NASA grant HST-AR-09927.01-A from the Space Tele- 
scope Science Institute, which is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS5-26555. 

REFERENCES 

Abramowitz M., Stegun I., 1974, Handbook of Mathematical Functions, Dover: New York 
Andredakis, Y.C., & Sanders, R.H. 1994, MNRAS, 267, 283 
Baes M., Dejonghe H., Buyle P., 2005, A&A, 432, 411 

Balcells M., Graham A.W., Dominguez-Palmero L., Peletier R., 2003, ApJ, 582, L79 
Begelman M.C., Blandford R.D., Rees M.J., 1980, Nature, 287, 307 
Bcrtin G., Ciotti L., Del Principe M., 2002, A&A, 386, 149 
Binney J., 1980, MNRAS, 190, 873 

Binncy J., Tremaine S., 1987, Galactic Dynamics, Princeton University Press: Princeton 
Caon N., Capaccioli M., Rampazzo R., 1990, A&AS, 86, 429 
Caon N., Capaccioli M., D'Onofrio M., 1993, MNRAS, 265, 1013 
Cardonc V.F., 2004, A&A, 415, 839 

Cardone V.F., Piedipalumbo E., Tortora C, 2005, MNRAS, 358, 1325 
Ciotti L., 1991, A&A, 249, 99 

Cox D.P., 2000, Allen's Astrophysical quantities, New York: AIP Press; Springer, p. 341 
Dehnen W., 1993, MNRAS, 265, 250 
de Jong, R.S. 1996, A&AS, 118, 557 

Demarco R., Magnard F., Durret F., Marquez, I., 2003, A&A, 407, 437 
de Vaucouleurs G., 1948, Annales d'Astrophysique, 11, 247 

de Vaucouleurs G., 1959, in Handbuch der Physik, ed. S.Fliigge, Springer, Berlin, p. 275 
de Vaucouleurs G., Capaccioli M., 1979, ApJS, 40, 669 

Erwin P., Graham A.W., Caon N., 2002, in "Carnegie Observatories Astrophysics Series, Vol. 1: Coevolution of Black Holes 

and Galaxies," ed. L.C. Ho (Pasadena: Carnegie Observatories) 
Faber S.M., et al., 1997, AJ, 114, 1771 
Ferrarese L., Merritt D., 2000, ApJ, 539, L9 

Ferrarese L., van den Bosch F.C., Ford H.C., Jaffe W., O'Connell R.W., 1994, AJ, 108, 1598 
Gebhardt K., et al., 1996, AJ, 112, 105 
Gebhardt K., et al., 2000, ApJ, 539, L13 

Goldstein, H., 1980, Classical Mechanics, Addison- Wesley, Boston 

Graham A.W., 2001, AJ, 121, 820 

Graham A.W., 2004, ApJ, 613, L33 

Graham A.W., Driver, S., 2005, PASA, 22(2), 118 



18 Terzic & Graham 



Graham A.W., Erwin P., Caon N, Trujillo I., 2001, ApJ, 563, Lll 

Graham A.W., Erwin P., Caon N, Trujillo I., 2003a, in Galaxy Evolution: Theory and Observations, RevMexAA (SC), eds., 
V. Avila-Reese, C. Firmani, C.S. Frenk, & C. Allen, vol.17, p. 196-197 

Graham A.W., Erwin P., Trujillo I., Ascnsio Ramos A., 2002, in "Carnegie Observatories Astrophysics Se- 
ries, Vol. 1: Cocvolution of Black Holes and Galaxies," ed. L.C. Ho (Pasadena: Carnegie Observatories: 
http: / /www. ociw.edu/ociw/symposia/scrics/symposiuml/procccdings. html I 

Graham A.W., Erwin P., Trujillo I., Asensio Ramos A., 2003b, AJ, 125, 2951 

Graham A.W., Guzman R., 2003, AJ, 125, 2936 

Graham A.W., Laucr T., Colless M.M., Postman M., 1996, ApJ, 465, 534 
Hansen S.H., Moore B., 2005, MNRAS, submitted lastro-ph/0411473t 
Hernquist L., 1990, ApJ, 356, 359 
Hjorth J., Madsen J., 1995, ApJ, 445, 55 
Jaffe W., 1983, MNRAS, 202, 995 

Kandrup H.E., Mahon M.E., 1994, Phys. Rev. E 49, 3735 

Kandrup H.E., Sideris I.V., Terzic B., Bohn C.L., 2003, ApJ, 597, 111 

Kawano Y., Oguri M., Matsubara T., Ikeuchi S., 2004, PASJ, 56, 253 

Kormcndy J., 1985, ApJ, 295, 73 

Kormendy J., Djorgovski S., 1989, ARA&A, 27, 235 

Landau, L.D., Lipshitz, E.M., 1976, Mechanics, Buttcrworth-Heinemann, Burlington, MA 
Laucr T.R., 1985, ApJ, 292, 104 
Lauer T.R., et al., 1995, AJ, 110, 2622 

Lima Ncto G.B., Gerbal D., Marquez I., 1999, MNRAS, 309, 481 
MacArthur, L.A., Courteau, S., & Holtzman, J. A. 2003, ApJ, 582, 689 
Magorrian J., et al., 1998, AJ, 115, 2285 
Makino J., Ebisuzaki T., 1996, ApJ, 465, 527 
Makino J., Funato Y., 2004, ApJ, 602, 93 

Marquez I., Lima Ncto G.B., Capelato H., Durrct F., Gerbal D., 2000, 353, 873 

Marquez I., Lima Ncto G.B., Capelato H., Durrct F., Lanzoni B., Gerbal D., 2001, A&A, 379, 767 

Mazure A., Capelato H.V., 2002, A&A, 383, 384 

McLurc R.J., Dunlop J.S., 2002, MNRAS, 331, 795 

Mellier Y., Mathez G., 1987, A&A, 175, 1 

Merritt D., 2004, in "Growing Black Holes: Accretion in Cosmological Context" jastro— p h/0409290 I 
Merritt D., Fridman T., 1996, ApJ, 460, 136 

Merritt D., Navarro J.F., Ludlow A., Jenkins A., 2005, ApJL, 624, L85 
Merritt D., Poon M.Y., 2004, ApJ, 606, 788 
Merritt D., Valluri M., 1996, ApJ, 471, 82 

Milosavljevic M., Merritt D., 2001, ApJ, 563, 34 

Milosavljevic M., Merritt D., 2005, in "Living Reviews in Relativity" jastro-ph/0410364) 
Milosavljevic, M., Merritt, D., Rest, A., & van den Bosch, F.C. 2002, MNRAS, 331, L51 
Navarro J.F., Frenk C.S., White S.D.M., 1996, ApJ, 462, 563 
Navarro J.F., Frenk C.S., White S.D.M., 1997, ApJ, 490, 493 

Navarro J.F., Hayashi E., Power C, Jenkins A.R., Frenk C.S., 2004, MNRAS, 349, 1039 
Poon M.Y., Merritt D., 2002, ApJ, 568, L89 
Poon M.Y., Merritt D., 2004, ApJ, 606, 774 
Prugnicl Ph., Simien F., 1997, A&A, 321, 111 

Ravindranath, S., Ho, L.C, & Filippenko, A.V. 2002, ApJ, 566, 801 
Rest A., ct al., 2001, AJ, 121, 2431 

Sersic J.-L., 1963, Boletin de la Asociacion Argentina de Astronomia, vol.6, p. 41 
Sersic J.L., 1968, Atlas de galaxias australcs 
Simonneau, E., Prada, F., 2004, RevMexAA, 40, 69 
Siopis C.V., Kandrup H.E., 2000, MNRAS, 319, 43 

Stiavelli M., Miller B.W., Fergus on H.C., Mack J., W hitmorc B.C., Lotz J.M., 2001, AJ, 121, 1385 
Terzic B., 2003, ApJ, submitted jastro-ph/ 0305005 I 
Terzic B., Kandrup H.E., 2004, MNRAS, 347, 957 

Tonry J., Blakeslee J.P., Ajhar E.A., Dressier A., 1997, ApJ, 475, 399 
Tremaine S., et al., 1994, AJ, 107, 634 
Trcnti M., Bertin G., 2005, A&A, 429, 161 

Trujillo I., Asensio Ramos A., Rubino-Martin J. A., Graham A.W., Agucrri J.A.L., Cepa J., Gutierrez CM., 2002, MNRAS, 
333, 510 

Trujillo I., Erwin P., Asensio Ramos A., Graham A.W. 2004, AJ, 127, 1917 
Young C.K., Currie M.J., 1994, MNRAS, 268, Lll 
Zhao H.S., 1996, MNRAS, 278, 488 



Density-potential pairs for stellar systems 19 

7 APPENDIX A: MASS, CIRCULAR VELOCITY, AND VELOCITY 
DISPERSION 

7.1 Prugniel— Simien density profile 

In this section we provide some additional helpful expressions related to the density profile 
p(r) given in Prugniel & Simien (1997, their equation B6), and given here by equation (@J). 
Assuming spherical symmetry, the enclosed mass, M(r), is simply 

r 

M(r) = Air J p(f)f 2 alf. (24) 
o 

Using the change of variable Z = b(r / Re) 1 /™ , such that df = R c nb~ r ' ''Z n ~ A dZ ', one obtains 
M(r) = 47rp #c 3 ^ n(p " 3) 7 (n(3 -p),Z), (25) 

where Z = b(r/R e ) l / n and 7(0,2;) is the incomplete gamma function given in equation ©. 
The total mass is obtained by replacing 7 (n(3 — p),Z) with T (n(3 — p)). 
The circular velocity is given by 



„ dIo(r ) = (26 ) 

Assuming isotropy, the spatial velocity dispersion is given in Binney (1980) as 

11 \ G f°° ,-, M ( f )j- 

°- (r) = W)L p{r) — dr 

p[r) Jz v ' 

This expression is easily computed numerically after making the substitution Z = [Zj cos 8), 
such that dZ/dO = Zsm.8/ cos 2 8, giving 

„ s 4irG(p R e nb n (P-V) 2 Z- n (P + V r /2 
a 2 r) = ^ e 



tan^(cos^) n(p+1) e- z/cos9 7 (n(3 -p), -) d9. (28) 

Jo V cos 8 J 



p(r) 

The projected, line-of-sight velocity dispersion as given by equation (B12) in Prugniel 
& Simien (1997) is 

o~p {R) = TfI3 ^, / p(f)M(f) df 



I(R)M tot J f 
87TGp 2 R c 3 n 2 b n ( 2p -V 



I{R)M X 



tot 



J™ ^Jz 2n - Z 2n Z- n ^- 1 eT* 1 (n(3 - p), z) dZ 



Jo v V cos 8 J 



I(R)M tot 



with Z now equal to b(R/R c ) 1/n . 

The mass profiles, circular velocity profiles, spatial and projected velocity dispersion 
profiles are shown in Fig. for a range of profile shapes. They agree closely with those 
obtained from the exact deprojection of Sersic's i? 1//n model (Ciotti 1991; Simonneau & 
Prada 2004). In order to show the behavior of f C j rc near the break radius r&, and to allow 
comparisons of (iw) 2 with a 2 , we present both t> circ on a logarithmic scale and (f C irc) 2 on a 
linear scale. 
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Figure 11. Core-less galaxies. Panel a) mass ('equation 1251 . b) circular velocity (equation 1261 . c) circular velocity squared, d) 
spatial velocity dispersion squared ('eouation l27l . and e) the square of the projected, line-of-sight velocity dispersion feauation l29l . 
for varying values of the profile shape n: n = 0.5 (solid lines), n = 1 (dotted), n = 2 (dashed), n = 4 (dash-dot), n = 10 (double- 
dash). We have set i? e =l, Mtot = 1, and G = 1. The square of each velocity term is thus effectively normalised to GMtot / Re- 
Spherical symmetry and isotropy have been assumed. 



7.2 Our density profile 

Below are the equivalent expressions relevant to our adaptation of Prugniel & Simien's 
density profile for galaxies with power-law cores. 
The enclosed mass M(r) is 

r r min r 

M(r) = An J p(f)f 2 df = An J pi(f)f 2 df + An J p 2 (f)f 2 df = Mi(r) + M 2 (r) (30) 

r min 

where pi(f) and p 2 (j) are given in equation (fTTj) . r m ; n is the minimum value of r or r?,, and 

M 1 (r)=AnJ 1 (r) J (31) 

M 2 (r) = 4vrL 2 (r), (32) 

where Ji(r) and L 2 {r) are defined in equations © and (|T2"j). respectively. The total mass 
is obtained by setting the second incomplete gamma function in the expression for L 2 (r) 
( equation 112)) equal to zero. Assuming spherical symmetry, the circular velocity is given by 
equation (j2HJ). 

Again assuming isotropy, the spatial velocity dispersion is given by 
cr s 2 (r) = K x {r) + K 2 (r) + K 3 (r), (33) 
where 

Ki{r) = — - / Pi{r)— —dr 

AnGp b W" \ 3^ 111 if ^ = 1 

Pi(r) \ 2(33^) (w 2 ' 1 "^ - r 2 M) if 7 ^ 1 and < 3 ' ^ 

K 2 (r) = -^-Mn) r ^-df 
p 2 (r) Jr max r 2 

= GMMab^ i~ z-^^-itf (35) 
R e p 2 {r) Jz max 

. ^ G r 00 . N M 2 (f) , 

K,{r) = —- P2(f)^df 

p 2 {r) Jr m „ r 

AnGp b 2 p 2 R 2 n 2 b 2n ^ 
P2{r) 
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x f°° Z^ p+1) ~ 1 e- 2 [T(n(3-p),Z b )-T{n(3-p),z)]dZ, (36) 

with r max = max(r,r 6 ), Z = b(r/R c ) 1/n , Z b = b(r b /R c ) 1/n and Z max = max(Z,Z b ). In 
passing, we note that the term — n(p + 1) appearing in the expression for K 2 (r) equals 
0.6097 — In — 0.05563/n, and is thus negative for values of n greater than about 0.2. For 
this reason we did not express K 2 {r) in terms of the incomplete gamma function. Following 
what was done with equation (J2*7|) . to help evaluate K 2 (r) and K${r) numerically one may 
apply the change of variable Z = Z max / cos 6, to give 

= GM^Pnb^ (p+1) r^^^VW-^ (37) 

R c p 2 {r) Jo 

_ ^G Pb YR e W^ {p+1) 
Aslrj - ^ ^ max j 

x r /2 tan^(cos^) n(p+1) e- Zmax/cos9 [r (n(3 —p),Z b ) —T (n(3 - p), Z max / cosflWflfi.) 

Alternatively, one may use the change of variable f = 1/s to transform the K 2 (r) term into 
the expression 

# 2 (r) = —-M^n) / p 2 (l/s)ds. (39) 



P2(r) 

The projected, line-of-sight velocity dispersion is found to be 
a v \R) = Si(R) + S 2 (R) + S 3 (R), (40) 
where 



Si(iJ) = r , m ,, / V 2 pi(f)M!(f) df 



I(R)M tot J f 5 

it 

%ixGoi?r h 2l R 3- 2 7 n 5n(2 7 -3) r ^ max ,- 

^£(3-,) L ^^Z^-HZ if 7 <3, (41) 



2G 



I(R)Mi 



tot 



Mi(r 6 )n6 n >p f° V^ 2 ™ - Z^Z-^+^-V^Z, (42) 



2(7 7 \/f 2 - i? 2 

*W - ITrSmZ I ^-^PMr) df 



I{R)M U 
8nGp b 2 p 2 R c 3 n 2 b n V p -V 



'tot 

J ' ] 1 1 ; i x 



I(R)Mi 



tot 



x f° ^Z 2 ™ - ^^-«(p+i)-i e -« [r (n(3 - p), Z 6 ) - T (n(3 -p), z)] dZ, (43) 

" ^raax 

with = r 6 , # max =m&x(R,R b ), Z = b(R/R c ) 1/n and Z fe = b(R b /R c fl n . Equation M and 
1321 can be integrated numerically using the change of variable Z = Z max / cos 6. 

The mass profiles, circular velocity profiles, spatial and projected velocity dispersion 
profiles are shown in Fig. ^] for a range of profile shapes. The intriguing "dip" in some of the 
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Figure 12. Core galaxies. Panels a-d) show the normalised mass (equation |2SJ , e— h) circular velocity feauation l26l . i— 1) circular 
velocity squared, m-p) spatial velocity dispersion squared fcquation l33l and q-t) the square of the projected, line-of-sight velocity 
dispersion ( equation 1401 . for varying values of the outer profile shape n: n = 0.5 (solid lines), n = 1 (dotted), n = 2 (dashed), 
n = 4 (dash-dot), n = 10 (double-dash); and varying central cusp slope 7. We have set R c — 1, r;, = 0.01, Mtot = 1, and G = 1. 
The square of each velocity term is thus effectively normalised to GM tot / R c . Spherical symmetry and isotropy have been assumed. 
To show the behavior of u c i rc near rj, and to enable easy comparison of (u c irc) 2 with and a^, we present v c - lIC on a logarithmic 
scale and (i> c i rc ) 2 on a linear scale. 

velocity dispersion profiles, most notably the high n profiles with small 7, is perhaps worth 
commenting on. It arises from the use of a sharp transition in the density profile (equation IHJ) 
when the inner and outer slope (on either side of r&) are markedly different. Application of 
a model with an extended transition region quenches this dip. The use of equation El with 
decreasing values of a, i.e. an increasingly broader transition region, steadily reduces and 
eliminates the dip. While profiles with a = 10 largely reproduce the velocity dispersions 
shown, density profiles with a = 1 have smooth velocity dispersion profiles. In any case, 
an inspection of the "pseudo" potential (equation [THJ) reveals that such a velocity structure 
does not give rise to an unstable situation for particle orbits. 

Equations pertaining to Prugniel & Simien's (1997) density profile are obtained by setting 
rfe = and pbf) = Pq. Indeed, upon making these substitutions in the derivation of equations 
for enclosed mass and velocity dispersion above, r min = 0, r max = r, i? max = R, Z mSuX = Z 
and the terms i^i(r), K 2 (r), Si(r) and S 2 (r) vanish. Equations (j2SJ), (J2*T|) and (J2~9"j) are then 
easily recovered. 
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8 APPENDIX B: DERIVATION OF THE POTENTIAL, FORCE AND 
INTENSITY PROFILES 

8.1 Potential 

If the expression for the density (equation EJ) is rewritten as p(r) = p\(r) + p 2 {r) ( a s was 
done in subsection 12.31 equation ITTj). then the potential (equation Ej) can be written as 



$0) 



-4vrG 



$(r) = -AttG 



r r;, oo 

/ Pi(r)f 2 df + J pi(f)f df + / p 2 (f)f df if r < r&, 

r r b 

r;, r oo 

- / pi(f)f 2 df + - / p 2 (f)f 2 df + J p 2 {f)f df if r > r b , 
o n r 

$Ji(r) + J 2 (r) + L 1 (r b ) if r < r 6 , 
±Ji(r 6 ) + ±L 2 (r) + Li r) if r > r 6 , ' 



(44) 



(45) 



The derivation of integrals J\{r) and J 2 {r), which only involve a power-law integrand as- 
sociated with p\(r), is trivial, while the integrals which involve p 2 (r) require a change of 
variable Z = b{f/R e ) 1 ' n , so that f = R c {Z/b) n and df = R e nZ n ^ n dZ: 



Loir) 



p 2 (r)r dr = p b p 
PbpRc 



Re 



1-p 



R r 



Pb 



pR e 2 nb n( - p ~ 2) 



p b pR c 2 nb n(p - 2) T (n{2-p),b 



f n(2- P )-l e -Z ^ 



^ ^ l/n N 

Re 



b(r/R a ) 1/r 



p 2 (r)r 2 dr = p b p 



r 

Rt 



-b{r/R e ) 1 ' n f 2 df 



n 



PbpRc 



i_y %-6(f/«.)v- £ 



fe(r/fl c ) 1/n 



Pb 



pR e 3 nb n( - p ~ 3) 



z n(3-p)-l e -Z dz 



p b pR e 3 nb n( - p - 3) 



r[n(3-p),b( r £ 



l/n N 



(46) 



r n(3-p),6 — 



l/n N 



(47) 



8.2 Force 

Computing the radial force requires differentiation of the potential (jHJ) with respect to the 
radial coordinate r. When r < r b , the differentiation is trivial after realising that the last 
term is constant. For r > r b , we use the relation given in Abramowitz & Stegun (1974, their 
equation 6.5.25): 

dT ^' x) = -x^e-*, (48) 
ox 
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and we note the following cancellation of terms simplifies matters 

lJ_Lp_ = J_LM_ = pbpRc (l-Y^e-Kr/^ 
r dr dr \R J 



(49) 



8.3 Projected Intensity 

Substituting pi(r) from equation (JT7J) into equation (fTBj) yields 



h{R) 



T 



Pbn' 



,1-7 



R 



Vr 2 - R 2 



dr. 



(50) 



By using the change of variable y = r 2 , it is trivial to compute the expression for the special 
case when 7 = 0. The solution is given in the first line of equation fi"9|) . 

When 7 = 1, the change of variable r = Rj cos 9, such that dr/d9 = Rsm9/ cos 2 9, gives 



h(R) 



r pbrb " 1 7^=W dr = r pbn " I ZoJ9 

R 



d9 



2 / 1 sin^ 

—pbn 1 In — r + — - 
I \ cos 9 cos y 



cos (R/ri,) 



2 _ /l + VI- cos 2 

^pbH 1 In - 

T \ cos 9 



cos 1 (R/ri > ) 



—pbTb 1 In 



1 - R 2 /r b 2 



R/n 



—pbn 1 In 



r & + y r b - R2 

R 



(51) 



For other values of the cusp slope 7, further transformation of equation ()50|) is required. 
Using the change of variable y = 1 — R 2 /r 2 , equation (jHUJ) becomes 



h(R) 



jPbTb 



pbn 



2 



l-R 2 /r b 



R 1 J B^ R 2 /rb 2 (-, 



(52) 



When 7 = 2, the incomplete beta function (equation 120)1 can be simplified through the 
change of variable u = t 2 , to give 



h(R) 



T 



pbn 



R- f * 
J 



dt 



2 -P^R-H\^[^l-R 2 /r 2 



(53) 



